Bayesian inference for high-dimensional linear regression under mnet priors

Abstract: For regression problems that involve many potential predictors, the Bayesian variable selection (BVS) method is a powerful tool, which associates each model with its posterior probabilities, and achieves superb prediction performance through Bayesian model averaging (BMA). Two challenges of using such models are, specifying a suitable prior, and computing posterior quantities for inference. We contribute to the literature of BVS modeling in the following aspects. First, we propose a new family of priors, called the mnet prior, which is indexed by a few hyperparameters that allow great flexibility in the prior density. The hyperparameters can also be treated as random, so that their values need not be tuned manually, but will instead adapt to the data. Simulation studies are used to demonstrate good prediction and variable selection performances of these models. Secondly, the analytical expression of the posterior distribution is unavailable for the BVS model under the mnet prior in general, as is the case for most BVS models. We develop an adaptive Markov chain Monte Carlo (MCMC) algorithm that facilitates the computation in high dimensional regression problems. Finally, we showcase various ways to do inference with BVS models, highlighting a new way to visualize the importance of each predictor, along with estimation of the coefficients and their uncertainties. These are demonstrated through the analysis of a breast cancer dataset. The Canadian Journal of Statistics xx: 1–??; 20?? c © 20?? Statistical Society of Canada Résumé: Insérer votre résumé ??? La revue canadienne de statistique xx: 1–??; 20?? c © 20?? Société statistique du Canada